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EXISTENCE OE WEAK SOLUTIONS TO A CLASS OF FOURTH ORDER PARTIAL 
DIFFERENTIAL EQUATIONS WITH WASSERSTEIN GRADIENT STRUCTURE 

DANIEL LOIBL, DANIEL MATTHES, AND JONATHAN ZINSL 


Abstract. We prove the global-in-time existence of nonnegative weak solutions to a class of fourth order 
partial differential equations on a convex bounded domain in arbitrary spatial dimensions. Our proof relies 
on the formal gradient flow structure of the equation with respect to the L^-Wasserstein distance on the space 
of probability measures. We construct a weak solution by approximation via the time-discrete minimizing 
movement scheme; necessary compactness estimates are derived by entropy-dissipation methods. Our theory 
essentially comprises the thin film and Derrida-Lebowitz-Speer-Spohn equations. 


1. Introduction 

1.1. Main results. This work is concerned with nonnegative solutions u : [0,oo) x SI —>■ R+ to the partial 
differential equation 

dtu{t, x) = V • u{t, x)V [Fz{x, u{t, x),Vu{t, x)) — V • Fp{x, u{t, x), Vu{t, x))] (1.1) 

for (t,x) € (0,oo) X SI, where SI C (d > 1) is a bounded and convex domain with smooth boundary dil 
and exterior unit normal vector field u. Additionally, the sought-for solution u is subject to the no-flux and 
homogeneous Neumann boundary conditions 

u{t,x)di, [Fz{x,u{t,x),Vu{t,x)) — V • Fp{x,u{t,x),yu{t,x))] = 0 = d^u{t,x) (1.2) 

for t > 0 and x £ 9S1, and to the initial condition 

m( 0, •) = Mo G .b^(Sl), uq > 0 and / Mo(x)dx = l. (1.3) 

Jn 

For the nonlinearity F - to which we refer as Lagrangian ~ we assume 
Assumption 1.1 (Conditions for F). F £ C^{Q x 1R+ x R'^), with: 

(i) Radial symmetry; There is some F : ft x R+ x R+ —?► R with drF(x, z,r) > 0 such that F(x,z,p) = 
F{x, z, IpI) for all p G R"^. 
faj Convexity; There exists 7 > 0 such that 

D(x,z.p)^(a;,2,l')[(C,C,7r),(^,C,7r)] > yjTTp 

for all (5, Ci'^) £ R"^ X R X R'^. and all (x,z,p) £ SI x R+ x R'^. 

(Hi) Bounds; There are constants C > c > 0 such that 

c|pp < F(x, z,p) < Cdpp -f 1 ) 

for every x £ SI, z £ R_|_, p gW^. Furthermore, there exists D > 0 such that 

\Fx{x,z,p)\, tr(F 3 , 3 ,(x,z,p)), z\F^{x,z,p)\, |Fp(x, z,p)|^ < D(|p|^-h 1). 

Our main result is concerned with the existence of solutions to (fTTIl-lfOl): 
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Theorem 1.2 (Existence of weak solutions). Assume that F satisfies Assumvtion \l.l\ and let a nonnegative 
initial datum uq € H^{n) with ||uo||li = 1 &e given. Then, there exists a nonnegative weak solution 
u : [0,oo) X n — >■ [0,oo] to (ll.ip - (ll.3ll in the following sense: For all ip £ C°°(n) with = 0 on dVL and 
all rj £ C))°((0, oo)) with 77(0) = 0, one has 


where 



u{t, x)(p{x)dtr]{t) dx dt 




X, u(t, x), <p(x))'g(t) dx dt, 


(1.4) 


N{x,p,il}) := F3,(x, p, Vp)- V-!/'+(E(x,/o, Vp)-pF^(x,/9, Vp))A'0 + Fp(x,/o, Vp)-(V/9A7/; + /9VA'0 + V^V'V/o). 

Partieularly, u £ L°“([0, T]; i?^(fl)) fl L^([0, T]; i?^(n)) for each T > 0, and u(t,-) is a probability density 
on fl at each t>0. Furthermore, the map t J^F{x,u{t,x),'Vu{t,x))dx is almost everywhere equal to a 
non-decreasing function. 


The second part of this paper is devoted to a special case of equation dm where the Lagrangian F does 
not satisfy the conditions from Assumption 11.11 Specifically, we let 

F(x,z,p) := i/'(z)2|pp, (1.5) 

for a map / : M+ —>■ K subject to 


Assumption 1.3 (Conditions for /). 

(i) f £ C^((0,oo)), /(O) = 0, f"{z) < 0 for all z > 0 and there exist C > 0 and a £ ([5 ~ ; l) such 

that for all z > 0.' f'{z) > Cz°‘~^. 

(ii) The limit \iTazf'{z) exists. 

z\0 

(Hi) There exists 5 > 0 such that ^ > 5 + 1 — ^ + + 8d for all z > 0. 

A typical example satisfying Assumption 1 1. 31 is the square root f{z) = s/z — see also Remark iTAl below. 
Note that since f is assumed to be unbounded at zero, especially the convexity assumption (jn| is not met 
by F. Nevertheless, we have 


Theorem 1.4 (Existence of weak solutions: non-convex case). Assume that F is as in (11.51) . and that f 
satisfies Assumvtion M.A Let moreover a nonnegative initial datum uq with ||mo||li = 1 /(^o) £ 

be given. Then, there exists a nonnegative weak solution u : [0,oo) x LI ^ [0, 00 ] to (11.11) - (11.31) ; For all 
ip £ C°^{Ll) with dvP = 0 on dfl and all rj £ (^((“((0, 00 )) with r]{0) = 0, one has 



u(t, x)p{x)dtr/{t) dx dt 



Aff (u(t, x), p{x))r](t) dx dt, 


where 

Affip, := A/(p)V/(p) • ViA + p/'(p)A/(p)AV;. 

Particularly, one has u £ L°“([0, T]; L^’(fl)) for some p > 1 and f{u) £ L°“([0, T]; iJ^(f2)) D A^([0, T]; F[‘^(Ll)) 
for each T > 0; and u{t,-) is a probability density on LI at each t > 0. Furthermore, the map t 1 —>■ 
fn almost everywhere equal to a non-decreasing function. 


1.2. Strategy of proof. The cornerstone of the proofs of Theorems II.21 and II.41 is the variational structure 
of equation (11.11) : Formally, it can be written as a gradient flow of the free energy functional $ : ^{Ll) —>■ 
M U {+ 00 } defined via 


$(p):= 


F[x, m(x), Vtt(x)) dx, 
+ 00 , 


\i p = u - and u £ Fd^lLl), 
otherwise. 


( 1 . 6 ) 


with respect to the L^-Wasserstein distance W 2 on the space of probability measures 3^{Ll) on LI. Indeed, 
for sufficiently regular F and u, one finds 


dtu = V • uV 


^<i>(tt) 

5u 
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where is the first variation of $ in at u. Since $ is finite on a subspace of absolutely continuous 
probability measures only, we shall from now on (by a slight abuse of notation) identify the measure /i with 
its Lebesgue density u and also write u G 

Note that the energy functional $ is in general not A-convex along geodesics in (,^(17), W 2 ) Q- Hence, 
the theory developed by Ambrosio, Gigli and Savare [l| is not immediately applicable. Still, a continuous 
flow for dni) can be defined as the time-continuous limit of the so-called minimizing movement scheme, 
dating back to De Giorgi 0: 

Fix a step size r > 0 and define a sequence (u")ngN recursively by 

:= uo, u” G argmin ( ^W 2 ('U, -I-$(«) ) (n G N). (1.7) 

\zT } 

With this sequence, we define the discrete solution Ut as the piecewise constant interpolant, i.e. for each 
t > 0, we set 


Ur{t,-) = u'^ forn=|’t/r]. (1.8) 

The strategy of proof for Theorems ll.2l and ll.4l can now be summarized as follows: First, we prove that the 
scheme (ll.7l) fc (ll.8ll is well-defined, that is, the successive minimizers in (11.711 exist. By an application of the 
flow interchange technique from [20j, we then derive an additional regularity property for the time-discrete 
solution Ur- In the framework of Assumption 11.11 the respective entropy-dissipation estimate formally 
amounts to 

— ^ f ulog(u)dx >7 f IIdx — C. (1.9) 

Jn Jn 

We consequently are in position to derive an approximate weak formulation satisfied by the time-discrete 
solution curves, and can pass to the continuous-time limit r 0 afterwards recovering the weak formulation 
of the original problem (jl.4|) . 


1.3. Examples and relation to the literature. The method used here is by now almost a standard 
technique to construct weak solutions to equations with a formal gradient flow structure. Since the seminal 
paper by Jordan, Kinderlehrer and Otto |15j | on the variational structure of the Fokker-Planck equation, 
various second order equations (e.g. [^1. fourth order equations (e.g. [H, [110,0) and systems (e.g. 
[30,0,01) have been analysed using this strategy. As models for physical or biological processes, sensible 
solutions should admit nonnegative values only. Since there exists no general comparison principle for 
equations of fourth order, this is a nontrivial issue. Thus, considering the dynamics on a suitable space of 
nonnegative measures, nonnegativity can be directly obtained from existence — and a formal Wasserstein 
gradient flow structure mi^t be an indication for it. However, in more specialized situations, other methods 

could be used (see e.g. [Si, 0,0,0). 

Our framework comprises two important examples: 

• The thin film equation 

dtu = —V • mV Am, (1-10) 

modelling the spreading of a liquid film over a solid surface, can be interpreted as the evolution 
of the Wasserstein gradient flow of the Dirichlet energy 4>(u) = ^||Vm ||^2 [20|, which fits into the 
framework of Assumption 1 1 . 1 1 for F{x,z,p) = Equation (11.101) can also be considered in the 

more general form 

dfU = —V • m(u)VAM, (1.11) 


with a possibly nonlinear mobility function m : R+ —>■ For the physically relevant cases, the 

existence of solutions to (11.111) has been shown e.g. in 0, IS]. Note that - with respect to a modified 
version of the Wasserstein distance - a formal eradient structure of (11.111) mav also be available ra¬ 
in the simplified framework of one spatial dimension, the long-time asymptotics of the classical thin 
film equation ()1.10|) have been analysed in 0], also with entropy-dissipation methods. 

• The Derrida-Lebowitz-Speer-Spohn equation 



dtU = V • uV 


( 1 . 12 ) 
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arises in the description of interface fluctuations in the Toom model 0[il as well as in models for 
semiconductors (see e.g. 16|). Often written in the form of the quantum drift diffusion equation 

( 


dtu = V • mV 


V 


V 


with a confinement pot ential V : ^ R, existence and long-time behaviour of solutions were 

studied e.g. in [1, EM Q • In QM. the formal Wasserstein gradient flow structure of (11.121) 
with respect to the Fisher information functional $(u) = 4|| Vi/m ||^2 is employed to construct weak 
solutions via the minimizing movement scheme. 

Actually, one may consider the more general class of functionals <1>/3 (m) = m^| Vm^ da; for /3 < 0, which are 

in general not geodesically A-convex with respect to the Wasserstein distance. Note, however, that in the 
range jd € [-5,-4], displacement convexity can be proved, see 0]. Theorem 11.41 presented here especially 
comprises the range /3 € [—1,0), by considering f{z) = z°‘ with a = 1 + §, which is in accordance to 
Assumption 11.31 Thus, this work provides a generalization of the existence results from for a wider class 
of energy functionals. 


1.4. Outline of the paper. After a brief preliminary section, we prove the well-posedness of the minimizing 
movement scheme under Assumption 11.11 in Section [31 After that, the additional regularity of the discrete 
solution is deduced in Section |4l and a discrete weak formulation is derived in Section |5l The proof of 
Theorem 11.21 is then completed by passing to the continuous-time limit (Section |6|). The extension of the 
strategy to the framwork of Assumption 11.31 will be sketched in Section [71 


2. Preliminaries 

In this section, we briefly summarize basic facts about gradient flows in the space of probability measures 
on the (bounded) set fl. For more details on optimal transport and gradient flows, we refer to the monographs 
by Villani [ 2 ^ and Ambrosio et al. [l|. 

A sequence {p,n)neN in .^(17) is said to converge narrowly to some limit probability measure /i G 
if for all continuous and bounded maps / : 17 —>■ M, one has 

lim f f{x)dfj,nix)= [ /(x)d/j(x). 

Note that if all measures are absolutely continuous and their respective densities converge weakly in LP{n) 
for some p > 1, narrow convergence of to p follows. 

The space ^^(17) can be endowed with the so-called L'^-Wasserstein distance W 2 defined as 

W 2 (/ro,/ii) := f inf [ \x - y\'^ djix^y)) 

\7er(/io.Mi) Jnxn / 

where r(/ro, pi) denotes the set of all transport plans from po to pi, i.e. 

r(/io, /ii) := {7 G X 17) : 7 has marginals po and pi} . 

Since 17 is a bounded subset of the infimum above is attained without further restrictions (e.g. finiteness 
of second moments) on the measures po and /ii. Moreover, narrow convergence and convergence with respect 
to W 2 are equivalent. Furthermore, if po is absolutely continuous with respect to the Lebesgue measure, 
the minimizer for W 2 is given in terms of the push-forward pi = for a measurable map T : 17 —>■ 17, 

and the Wasserstein distance reads 

W 2 (/ro,/Ji)^= [ \x-T{x)\^ dp,Q{x). (2.1) 

Jn 

In this non-Euclidean framework, one can dehne a notion of convexity via the so-called displacement inter¬ 
polation: A functional 21 : J^(17) —>■ R U {+ 00 } is called displacement X-convex for some A G R if for each 
POiPi S .^(17), there exists a constant-speed geodesic curve /is : [0,1] —>■ J^(17) connecting po and pi such 
that 

2l(/is) < (1 - s)2l(/io) + s2l(/ii) - -^ 5(1 - s)W 2 (/io, pif. 

For displacement convex functionals, we shall use the following notion of gradient flow: 
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Definition 2.1 (k-Aows fl)- Let 21 : ^ RU{+c»} he proper, lower semicontinuous and k- displacement 

convex w.r.t. W 2 for some it G K. ^ continuous semigroup on {Vl),'W 2 ) satisfying the evolution 
variational estimate 

w) + ^Wl{5f{w),w) + 2l(S?M) < QL{w) 

for arbitrary w,w in the domain of and for all s>0, as well as the monotonicity condition 

2l(Sf M) < 2t(Sf H) VO < s < t 
for all w G is called K-flow or gradient flow of%. 


The cornerstone of the rigorous derivation of (HH) is 

Theorem 2.2 (Flow interchange lemma (^ . Thin. 3.2]). Let % he a proper, lower semicontinuous and 
displacement X-convex functional on W 2 ) and assume that there exists a X-flow S^. Let furthermore 

dl be another proper, lower semicontinuous functional on [VC) ,'W 2 ) suchthatY)ouY[^) C Dom(2l). Assume 

that, for arbitrary r > 0 and w G the functional ^W 2 (-, possesses a minimizer w on 

Then, the following holds: 

Qi{w) + tD^'^[w) + ^W2(w, w) < Qi{w). 

There, denotes the dissipation of the functional $ along the X-flow of the functional 21, i.e. 


T)^'i>[w) := limsup 

s\,0 


T(ic) — 'it[Sf[w)) 
s 


( 2 . 2 ) 


The following theorem provides an extension of the Aubin-Lions compactness lemma to our metric setting: 


Theorem 2.3 (Extension of the Aubin-Lions lemma [2J, Thm. 2]). Let IX be a Banach space and A : X —>■ 
[0, 00 ] be lower semicontinuous and have relatively compact sublevels in X. Let furthermore g : XxX —> [0, 00 ] 
be lower semicontinuous and such that s[u, u) = 0 for u,u G Dom(A) implies u = u. 

Let {Uk)k^n be a sequence of measurable functions Uk ■ (0,T) —?► X. If 

pT 


sup / 
feeN Jo 


A[Uk{t)) dt < 00 , 


rT-h 


lim sup / 
keN Jo 


s{Uk{t -b h), Uk[t)) At = 0, 
then there exists a subsequence that converges in measure w.r.t. t G (0,T) to a limit U : (0,T) —>■ X. 


(2.3) 

(2.4) 


3. The minimizing movement scheme 

This section is devoted to the study of the energy functional $ defined in (11.61) and of the associated 
minimizing movement scheme (ini). Throughout, we assume that uo G D H^(Vl) and that F satisfies 

the conditions from Assumption 11.11 

Proposition 3.1 (Properties of the energy). The energy functional <I> is proper, nonnegative and lower 
semicontinuous in (b^(n), W 2 ). Moreover, there exist Co > 0 and Ci > 0 such that 

«>(m) > Collujl^i-Cl WuG A^[n)nH^[n). (3.1) 

Proof. Obviously, Assumption 11.1 [Hull yields that 0 < $(u) < oo for all u G CA{Vl) (~l H^{Vl). For those u, 
Assumption ll.ltlml) gives c||Vm ||^2 < $(u). Using ||u||li = 1 and Poincare’s inequality 

\\u-C\Vl)-^f^,<K\\yur,\\l., 

we deduce (EH). ft remains to prove that $ is lower semicontinuous. Let W2(^n,M) —t 0 as n —>■ oo for 
a sequence [pn)n&i and a limit p, in CA[Vl). If liminf <f>(/x„) = -|-oo, there is nothing to prove. We thus 

n^cx) 

can assume without loss of generality that (‘b(/i„))„gN is bounded. Hence pn is absolutely continuous for 
every n £ N, and by Alaoglu’s theorem, there exists a subsequence (non-relabelled) on which ^ u in 
H^{Vl). Since F is jointly convex in [x,z,p) by Assumptio n 1 1.1 ifnll and bounded w.r.t. a; £ U, <i> is lower 
semicontinuous with respect to weak convergence in H^[Vl) |23l . Thm. 10.16], completing the proof. □ 
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The properties from Proposition 13.11 allow us to show the well-posedness of the minimizing movement 
scheme uni). For the sake of presentation, we introduce the notation 

d>,(.;u) := + 

for the Yosida penalization associated to for given r > 0 and v G 

Proposition 3.2 (Minimizing movement scheme), //r > 0 and v G with $(?;) < oo, then there exists 

a minimizer u G o/<&,-(•; u) satisfying <I>(u) < oo. 

Proof. Since also (s v) is nonnegative, infimizing sequences for this functional are bounded. The existence 
of a minimizer is now directly obtained from the properties in Proi30sition l3.il recall that Un ^ u narrowly 
implies in particular that W 2 (M„,'y) —W 2 (m, u), since C is bounded. □ 

From the scheme ini), the following classical estimates directly follow (see for instance Q): 

Proposition 3.3 (Classical estimates). Let r > 0, uq G fl and define (u”)„gN and Ur by 

(Hilie/lLHi). Then the following holds: 

$«) < d>(u™) < $(mo) Vn > to > 0, (3.2) 

^W2«-\02<2t<1>(uo), (3.3) 

n>l 

W 2 {Ur{t),Ur{s)) < \/2^{uQ){\t - s| + t) Vs, t > 0. (3.4) 


Later, we need the following classical result on boundary values: 


Lemma 3.4 (Boundary terms 
Vv ■ V^vv <0 on dLl. 


l4 Lemma 5.2]). Let v G with d^v = 0 on dfl. 


Then, one has 


4. Additional regularity of the discrete solution 

The main result of this section is concerned with the additional regularity of the minimizers in (II.7I) and 
reads 

Proposition 4.1 (Additional regularity). Let r > 0 and v G fl Then, a minimizer u G 

^(il) (1 (Tl) o/4>t-(-;?;) admits the estimate 


[ ||V^M|pda;<— [ {v\og{v)-u\og{u)) dx + C3{\\u\\HI +1) 
Jn ir Jq 

with 7 > 0 from Assumption^TTT^ and some constant C 3 > 0. In particular, u G 


(4.1) 


Proof. The main idea of proof is to apply the flow interchange lemma ()Theorem l2.2|) to the free energy 4>, 
with the following choice of the auxiliary functional: It is well-known (see for instance [Ei) that Boltzmann’s 
entropy £(v) = n log(i;) da; is a lower semicontinuous displacement 0-convex functional on (.^3^(12), W 2 ), 
and its corresponding gradient flow semigroup (see Definition 12.1[) is a solution to the heat equation with 
Neumann boundary conditions 


a.Sf(i;) = ASf(u) 

5j.Sf (i;) = 0 


on (0,oo) X D, 
on (0,oo) X dn. 


(4.2) 

(4.3) 


For the dissipation of $ at u along the smooth flow S^, we obtain (write Us := Sf(u) for brevity) 

--^$( ms ) = - j T>(x,z,p)F{x,Us,Vus)[{0,Aus,A\7us)]dx 

d „ 

— ^ ^ 5 \^Xi (O5 ^Xi '^s i ^Xi ^^s)] dx. 

Jn 


2=1 ' 
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Integration by parts yields with the canonical unit vectors ei,..., € 

d 


= X! / ^ p)F{x,Us,yue)[{ei,dj;,Us,dj;yus), {0,dj;,Us,dj;yus)] dx 

as i=i "'a 

d p 

^ ^ tig , VUs) [(0, ^^s)] ^2 dfj. 

i=i 

For the boundary term, we use the boundary condition (14.31) and Assumption ll.llH) ) to get 

d p 

^ ^ ^1 ^Xi^si Us)\^i 

i^i 

= - Fz{x,Us,Vus)Vus ■ V da - / F'p(ai, Us, Vits) • dtr 

J dCl J 50 

= -/ F;(x,Us, |Vus|)—^ • V^Usudcr. 

JdQ |Vu^| 

The last expression above is nonnegative since F^ > 0 thanks to Assumption 1 1.1 tUl) and due to Lemma 13.41 
By elementary linear algebra, we have that 

d p 

/ Df^,.p-)F{x,Us,VUs)[iei,da:iUs,do,,VUs),iO,do,iUs,do,iVUs)]dx 
i=l "'a 

d p 

— ^ ^ j [( 2 7 ^Xi Us , VWg ) , ( 2 Gj, Ox^ Us , ^ )] 

d p 

“X! / D(,,,^,p)A(a;,Ms,Vus)[(5e*,0,0),(ie,,0,0)]da; 

i=l 

f ■ ■ 1 

>1 '^Idxyusf dx - -\\ti Fz;x{x,Us,Vus)\\l'^, 

using Assumption 11.1 [Hil) in the last step. With Assumption 11.1 Ipnl) . one moreover has 

\tTFxx{x,Us,VUs)\ < D{\VUs\'^ + 1), 
so all in all there exists C > 0 such that 

-^^'(ms)>7 f llVVfdx-Cdlit^ll^i + 1) > 7 [ llV^udl^dx - Cdlitll^i + 1), 

ds Jq Jq 

where the last estimate holds since ||us||jji is nonincreasing w.r.t. s > 0 due to the properties of the heat 
flow and Lemma 13.41 Passing to the limit inferior as s 0 yields by lower semicontinuity that 


D^<i>(u)>7 / llV^uf dT-Cdlull^x + l). 
Jn 


(4.4) 

□ 


Combining (14.41) with the flow interchange estimate (12.21) obviously yields (14.1|) . 

5. The discrete weak formulation 

We are now in position to derive a discrete version of the weak formulation dEl: 

Lemma 5.1 (Discrete weak formulation). Let r > 0 and define the discrete solution Ur by (11.71) 0 (11.81) . 
Then, for all ip G with = 0 on dfl and all rj G C^((0, oo)) D C'(M+), the following discrete weak 

formulation holds: 

f [ Ur{t,x)p{x) ^'^^ ^ ^ ^ dx dt — f f Af{x,UT{t,x),p{x))T]r{t)dxdt 

lo Jn X Jo Jn 

< T||¥’llcdl7llco^(Mo), 

where J\f{x, •, •) is defined as in Theorem \1.2\ and r]T{s) := rj ([I-] r) for s > 0. 


(5.1) 
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Proof. We use the so-called JKO method [I^ (see also ES [ 2 ^) introducing a suitable perturbation of the 
successive minimizers in the scheme uni). Specifically, for fixed ip G C°°{n) with d^ip = 0 on dft, we denote 
by X(.-) : R+ x fl fi the smooth flow associated to the ordinary differential equation 

y(s) = \7<p{y{s)), for s > 0. 

Using the minimality property of u”, one has for all s > 0: 


= i mXs#K) - $«)) + ^ {wUXs#K, <-i) - W2«, <-!)). 


(5.2) 


s . .... 2rs 

The discrete weak formulation m is obtained by passing to the limit as s \ 0 in (15.211 . Since u" and u" ^ 
are probability densities, there exists an optimal transport map T : such that = T^uf. (see 

El), so consequently 

Wi(x,#<,<-1)-W2«,<-1)< [ (lx.-rp-iid-TpXdx 

JQ 

< [ {Xs{x) - x) ■ {Xs[x) + x- 2T{x))u'f dx, 

Jn 

using the elementary identity |ap — = {a — b) ■ {a + b) for a,b G in the last step. By dominated 

convergence and the definition of Xg, one has 

-W^«, <-!))=-- / X(p{x) ■ {T{x) - x)u'fdx. 
s\o Its t Jq 

Using the Taylor expansion ip{T{x)) — (fi{x) = V(p{x) ■ (T(x) — x) -b i(T(x) — x) • V^ip(x){T{x) — x) for some 
intermediate value x, the right-hand side above can be recast to 


[ V(/j(x) • (r(x)-xXdx < i / «-< ^)ipdx + ^\\p\\c2Wl{ 
T Jn r 2 t 


2 /^.n „.n —1 


)■ 


To pass to the limit in the first term of the r.h.s. in (15.21) . we introduce the volume distortio n 14 := 
det(VXs) > 0 associated to the vector field Xg, for small s > 0. The following identities hold (see [20|, §2] 
for a proof): 


_d 

ds 


U, = A^, 


s =0 


_d 

ds 


s =0 L 




(Xg) 


-1 


oXg 


= —XuAif — uX Aip — 


(5.3) 


Using the change of variables x = Xg{y), we can rewrite ^{Xgp,u^) as follows: 


$(A,#<)= / U(A„^,V(-^o(X,)-i)oA,)U,d2/. 


By dominated convergence, one consequently has 

lim i(<l>(X,#<)-<!>«)) 

s \,0 S 


^{x,z,p)Fix,u'f,Xu^) 

= [ A/'(x,<,(/7)dx, 

Jn 


ds 


s=0 






(^s) 


-1 


oX 


+ F(x,<,V<). 


ds 


s =0 


f". dy 


by applying (15.31) . 

Putting together, we have proved that 

0< [ A/'(x,<,(/5)dx-b - [ « - <“^)(/5dx-b ;^||v7||c2W2«,u”"^). 

Jn T Jn 

Since J\f{x, p) is linear with respect to p, repeating the calculations above for —p in place of p yields the 
converse inequality and hence 

1 
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We now introduce a nonnegative temporal test function 77 G oo))nC(]R+), multiply (15.41) with r-qinr) 

and sum over n G N to obtain 


-rM\c4v\\co^{uo) <^ [ + ^ [ Tr]{nT)M{x,u'^,ip)dx 


(5.5) 




using the total square distance estimate (1531) to simplify the last term. Rearranging above as 


[ r]{nT){u^-u^ ^)ifdx 



r]{nT) — ri{{n + l)r) 
r 


u^(pdx, 


and introducing r]T{s) ■= ^([ 7 ] t): the discrete weak formulation (15.11) is obtained by rewriting (15.5|) in 
spatio-temporal integral form, recalling the definition of the discrete solution Ur- For sign-changing test 
functions 77 , we decompose into positive and negative part and subtract the respective estimates (15.51) also 
to arrive at dnn). □ 


6. Passage to continuous time 


In order to pass to the continuous-time limit t 0, the following a priori estimates are required: 

Proposition 6.1 (A priori estimates). For each fixed T > 0, there exist constants Ci,C 2 > 0 such that for 
all T G (0,r), the discrete solution Ur defined via (11.71) 61 (11.81) admits the estimates 

l|Wr||L°°([0,T];ffi) < c*! O.'nd ||uT-||L 2 ([ 0 ,T];il 2 ) < (72. 

Proof. Since d>(uo) is finite, the first estimate is an immediate consequence of the energy estimate (13.21) from 
Proposition 13.31 and (13.11) from Proposition 13.11 For the second estimate, we define N := G N and use 
the additional regularity estimate ()4.1|1 from Proposition 14.II 


r\\UrrH.dt= r\\Ur\\hdt+ \\V\r\\h dt 
Jo Jo Jo 

N 

< T’lhr||i=»([o.T];/IQ + E ^ 


n—1 

N 




-(6«-^) - 6(0) + CardKr^i + 1) 


=iL7 

Using the elementary estimate — 7 < -zlog( 2 :) < z'^ for 2 ; > 0 in combination with the uniform estimate on 
Ur in L°°([0, T]; iJ^(U)), one obtains (for some constant (7 > 0) 

~ 1 

lk.||L^([O.T];iI=) < C + -lluolli^ + 

7 76 

which is a finite value, thanks to the assumptions on 12 and uq. □ 


We are now able to identify a candidate for the limit curve: 

Proposition 6.2 (Convergence). Let T > 0, a vanishing sequence of step sizes 0 (A: —^ 00 ) and the 

associated family of discrete solutions {urf.)k£n be given. There exists a (non-relabelled) subsequence and a 
limit map u G (7i/2([0, T]; (.^(12), W 2 )) (7 L°°([0, T]; iJi(O)) n ^^([o, T]; ij2(0)) such that for k ^ 00 : 

(a) Ur^ —>■ u uniformly w.r.t. t G [0,T] in (.^(12), W 2 ). 

(b) Mt-j, —>■ u m L^([0, T]; i7^(f2)). 

(c) Ur^if,-) u{t,-) in 77^(12) for almost every t G [0,r], as well as Ur^ —t u and ’Vur^, Vit almost 
everywhere in [0,T] x 12. 

(d) Urk u in U^([0, T]; 27^(12)). 

Proof. Part (jaj) is an immediate consequence of the Holder-type estimate (13.41) from Proposition 13.31 and a 
refined version of the Arzela-Ascoli theorem [l|. Prop. 3.3.1]. The weak convergence (0 is obtained from 
Alaoglu’s theorem and Proposition 16.II It remains to prove the strong convergence (0 which yields (jcj) by 






10 


DANIEL LOIBL, DANIEL MATTHES, AND JONATHAN ZINSL 


extraction of suitable subsequences. For the proof of (|b|, we apply Theorem 12.31 to the family Let 

X := and define yt: X —[0, oo] with 


Aip) := 



+ 00 , 


otherwise, 


which has relatively compact sublevels in X by the Rellich-Kondrachov compactness theorem. Furthermore, 
define g ; X x X —>■ [0, oo] by 


, „ I W 2 (p, if p, p are probability densities on fl, 

^‘->’•^'■=[+ 00 . otherwise, 


which is an admissible choice for the application of Theorem 12.31 Using the technique from [26| , one ver¬ 
ifies that hypothesis (EH) holds, thanks to the estimates from Proposition 13.31 Since hypothesis (ESD 
coincides with the second part of Proposition 16.11 we conclude with Theorem 12.31 that (on a subsequence) 

I|w||l 2 ([o,t];H 1 ) as fc —>■ oo by dominated 


converges to in in measure with respect to t G {0,T). Using the uniform estimate in 

L°°([0, Tj; from Proposition 16.11 one has \\urJ\L'^(io,T]im) 


convergence. With weak convergence and the Radon-Riesz theorem, we consequently have (0. 


□ 


To complete the proof of Theorem 11.21 it remains to verify that u is a weak solution to (10. 

Proposition 6.3 (Continuous-time limit). Let {Tk)keN be the subsequence on which the convergence prop¬ 
erties from Proposition 16.31 hold. Then, 


lim 

k—^cio 




(t + Tk)- Pr, 


Tk 


— u{t, x)dtrj{f) ) p{x) dx dt = 0, 


lim / / {M{x,Ur,,{t,x),Lp{x))prk{t)~J^{x,u{t,x),ip{x))rj{t)\dxdt = Q. 

Jq\ ) 


Moreover, for almost every t S [0,T], one has 


lim $(M^,(t,-)) = «>(M(t,-))- 

k—yoo 


( 6 . 1 ) 

( 6 . 2 ) 

(6.3) 


Proof. Since + Tk) — Vrk) converges to dtp uniformly, (16.11) is an immediate consequence of the 

convergence properties stated in Proposition [10 For the proof of (16.21) . we first note that 

|A/'(x,p,(p)| < C||(p||c2(l + |Vpp) Va; G ft, 

for some constant C > 0, thanks to the bounds on F from Assumption 11.1 tHI)) fc (Hull . By Proposition |10 we 
ha.ve J\f{x,Urk{t,x),(p{x))prk{t) N{x,u{t,x),ip{x))p{f) pointwise a.e. on [0,T] x fl. Convergence of the 
integral is obtained by Vitali’s convergence theorem as follows: Let e > 0. Since Urk —>■ u in L^([0, Tj; iL^(fl)), 
there exists ^ G ^0, 

G N, one has 


2C||vp|lc2||r/|lcO 


^ such that for all Borel sets Ac [0,T] x with £‘^+^(A) < 6 and all 


' da; dt < 


2 ^ 117^1102 11^1100 


Consequently, for those A and all fc G N, 

\M{x,Urk{t,x),'p{x))prk{t)\ dt < C||</j||c2| 


loo 


< e. 


2C||va||c2||p||oo, 

which allows us to apply Vitali’s convergence theorem to the family {J\f{-,Urk, T)VTk )fcgN- The proof of (16.31) 
can be obtained by a similar argument: Starting from the first part of (j^ from Proposition 16.21 we deduce 
that for every subsequence of {Tk)kenj there exists a sub-subsequence on which ^{urj_) converges to 

$(u) as k ^ oo. For the convergence, we employ Vitali’s theorem in almost the same way as above; the 
necessary uniform integrability is a consequence of Assumption 11.1 lilil) . □ 


Obviously, (I6.1l) fc (l6.2l) imply that u satisfies the weak formulation (11.41) in Theorem 11.21 The remaining 
assertion on the monotonicity of ^{u{t, •)) w.r.t. t is an immediate consequence of the discrete estimate (13.21) 
and the convergence (16.31) . 
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7. Extension to non-convex Lagrangians 

In this section, we sketch the possible extension of the strategy described above to the case covered by 
Assumption 11.31 For the sake of brevity, we skip most of the technical details. Recall that the free energy is 
defined as 

^|V/(M)pdx, 

if u G and /(u) G and ^(m) = +oo otherwise. 

We first summarize several properties on the nonlinearity / which are elementary consequences of As¬ 
sumption 11.31 

Lemma 7.1 (Properties of /). The following statements hold: 

(a) With C and a from Assumvtion \1.3\. one has f{z) > —z°‘ for all z >0. 

(b) There exists a constant Cq > 0 such that f{z) < Co{z + 1). 

(c) There exists a constant Ci > 0 such that zf'{z) < Ci{f{z) + 1) for all z > 0. 

We are now in position to construct a discrete solution via the minimizing movement scheme (ini): 

Proposition 7.2 (Minimizing movement scheme). //r > 0 and v G with f{v) G then there 

exists a minimizer u G o/<I>^(-;z;) satisfying f(u) G 

Proof. Let an infimizing sequence {uk)keN for $,-(•, i;) be given. Since $ is nonnegative, ^riuk', v) is bounded; 
hence 

W 2 {uk,v)<C and ^{uk) < C 

for some C > 0. Using Poincare’s inequality together with Lemma ITA] yields \\f{uk)\\H^ < C since 

£(0)-! f f{uk)dx<CoC{n)-^ [ (Mfe + l)dx = Co(/:(0)-i + l). 

Jq Jq 

Using Alaoglu’s theorem in combination with the Rellich-Kondrachov compactness theorem, one has (on a 
non-relabelled subsequence) that Wk '■= f{uk) converges to some w weakly in as well as strongly in 

L^{Tl) and pointwise almost everywhere on 17. By Prokhorov’s theorem, Ufc —t m in (l3^(17), W 2 ), on a further 
subsequence. Clearly, u = f~^{w), so f{uk) f{u) in iJ^(17). By weak lower semicontinuity of <&,-(•;t) 
(following by weak lower semicontinuity of ^||Vw||^ 2 ), it is easy to conclude that u is a minimizer. □ 

As before, the minimizer gains in regularity: 

Proposition 7.3 (Additional regularity). Let t > 0 and v G with f{y) G Then, a minimizer 

u G 1^(17) o/d>r(’;^^) admits the estimate 

/ l|V^/(u)||^da; < ^ [ {v\og{v)-u\og{u))dx, (7.1) 

Jq ot Jq 

where 5 > d is a constant depending only on the structure of f and the spatial dimension d. 

Proof. We proceed as in the proof of Proposition 14.11 Recall the definition of the heat flow from 
(IT^fcdg. Clearly, estimate (ED is obtained via the flow interchange lemma fTheorem l2.2l) if the dissipa¬ 
tion along the heat flow admits the estimate 

D^$(u) > (5 / ||V2/(u)||2dx. (7.2) 

Jq 

Denote, for brevity, fg := /(Sf (it)) (and analogously /' etc.) and Ug := Sf (it), for s > 0. Since Ug is smooth 
and strictly positive, one obtains 

S. ■.= -4-Hus)= f Afgf'Augdx= [ {Afgfdx- f AfgfJ\Vug\^ 
els JQ Jn Jn 

using the product rule Afg = /"jVitsp -f f'gAug. Since A/^V/s • u = 0 and V/s • < 0 on dQ (recall 

Lemma 13.411 . we have with integration by parts that 

^s-S [ ||V2/gf dx> (1-J) f WV^fgfdx- f V/g-V(i7(/g)|V/gn dx, 

J Cl J Cl J Cl 
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where <5 G (0,1) is a parameter to be specified below and i?(/s) := —-jjfyr- Note that \7H{fs) = L{fs)Vfs 

f'" f' -|-2f f 

with L{fs) = — “ Ifiyi “ ■ By Gaufi’ theorem, we first obtain 

0= [ V-(i/(/«)|V/,pV/,)dx= [ [iJ(/,)A/,|V/,p + 2iJ(/,)V/,-VV«V/,+L(/,)|V/,|4] dx, 

JQ JQ 

d . 


and consequently, for y := 


^s-S f ||VV«f dx>(l-5) f WV^fsfdx 

JQ JQ 


+ [ [H{fs){{l-x)Afs\Vfsf-2xyfs-y^fsyfs)-xL{fs)\yfs\^] dx. 
Jn 


(7.3) 


Defining R := V^/s — (which is traceless and symmetric) and choosing the Frobenius norm as matrix 


norm, we observe that 


iVVsf = tr(VV.V2/,) = + ||i7f. 


Insert this into (17.31) to find 

^s-S f liVV.f dx> f 

jQ JQ, 


1-6 . 


(A/«)2+ 1- 1 + - X i7(/,)A/,|V/. 


+ (1 - 5)||i?f - 2xH{fs)Wfs ■ RWfs - xL{fs)\Vfs 


(7.4) 


dx. 


Applying Young’s inequality to the first two terms on the r.h.s. of (17.41) and using Lemma lA.ll from Appendix 
El below to estimate the third and fourth term, we end up with 


^s-6 [ \\y^fsfdx>x [ H{fsf\yfs 

JQ JQ 


LUs) 

H{fs? 1 - <5 




2 2 


dx. (7.5) 


We conclude the proof by showing that there exists 6 € (0,1) such that the expression in square brackets in 
dLSl) is nonnegative. Indeed, since 


L{fs) fj'f's 




by Assumption 11.31 one has 

L{fs) 1 


H{fsf 1 - 


— 1 — — + - \/d^ + > -— 


2 2 


l-<5 


d 1 


(5 — (5((5—1 — — + -\/d?+d>d 


which is nonnegative choosing 6 sufficiently small. The desired estimate dUl) now follows by passage to 
s 0 and lower semicontinuity. □ 


Remark 7.4 (Power functions). In consistence with the results from ]2(\1 . our theory eomprises the important 
ease that f{z) = Cz'^ is a genuine eoncave power function, where the exponent a is allowed inside the range 


4 4 

Note that, in view of Assumvtion \773l one has 

1 1 

W1 -I- - > 

2 4 4 


3 1 

— — + 


(7.6) 


^ r(z)/'(z) 2-0 ^ _ 

d-2 d nzr 


hence Assumvtion \1.3\ is equivalent to (17.61) for this special choice of f. 

With all these prerequisites at hand, we derive a discrete weak formulation and a priori estimates satisfied 
by the discrete solution Ur- 
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Proposition 7.5 (Discrete weak formulation and a priori estimates). Let r > 0 and define the discrete 
solution Ur by (11.71) & (11.81) . Then, for all fi > 0, ip G C°°{Ll) with = 0 on dLl and all rj G oo)) (~) 

C(R+), the following discrete weak formulation holds: 

~ II^IIc'’^t<1’(mo) +/3 [ f - \v\t{ + — 1 dt 

Jo Jq t 

< f f f [ r?rA//(Mr,V^)da;dt (7.7) 

Jo Jq ^ Jo Jn 

< \\r]\\coKT^{uo) — P [ [ ^ — lii^(t^x)logUr{t,x)dxdt, 

Jo Jn 

where Mf is defined as in Theorem \1.4\ and k > d is such that —k 1 < V'^p(x) < kI for all x G LI. For 
h : IR.+ —>■ R, we denote hrfs) = /i(|"|']r). 

Furthermore, the following a priori estimates hold for each fixed T > 0: 

ll/(uT)||L“([0.r];ffi) < Cl, \\f {Ur)\\L^{[o,T]-,H^) < <^2 and ||nT||L“([0,T];LP) < C's, (7.8) 

for some constants Cj > 0 and some p > 1. 

Proof. The a priori estimates on Ur are obtained in the same way as in the proof of Proposition 16.11 Note 
that, thanks to Assumption 11.31 and Sobolev’s inequality, we have for all t G [0,T]: 


(O 

VI 

ft. 

for all p G ( 1, Of— 

V d 

\\Ur{t,-)\\LP < C 

for all p G (1, oo) 

\\Ur{t,-)\\LP<C 

for all p G (1, oo] 


2d 


d-2 


if d > 3, 


if d = 2, 
if d= 1. 


(7.9) 


To derive the discrete weak formulation dm), we use the flow interchange lemma iTheorem 12.21) with the 
displacement (—it)-convex regularized potential energy 


V{v) 


V log V dx 


vpdx, 


Jn Jn 

for P > 0. The associated (—K)-flow is given by (see [ij for more details) 

dsS^(v) = /3ASy (n) + V • (S^ (n)Vtp) on (0, oo) x Lt, 

di,S'^{v) =0 on (0,oo) x dLl. 


Let T > 0 and n G N. Since (u”) is smooth and strictly positive, we obtain for s > 0 (writing again 
Ws := Sy(ny) and /« := f{us) for brevity): 

-^d>(us) = [ /'(ms)A/(ms)(/3Ams + V • (itsVfp)) dx. 
ds Jq 

The viscosity term can be treated exactly as in the proof of Proposition 17.31 yielding 

--^^'(ms) = J [/3d||V^/s||^ + A/(us)V/s • V(/? + A/s/'(us)msA(/?] dx. 

Using Holder’s inequality and Lemma l7.1f ci. we can estimate from below: 

-^$(«.) > mi^^fsWh - CWV^fsh^iWfsWm + 1), (7.10) 

for some constant C > 0 depending on p. The following Sobolev inequality holds for all p G 

\\p\\m<C'\\V^p\\Up\\]j', with d = ^G (0,1). (7.11) 

Using (I7.11|) in (17.101) . recalling that |1 /s||li < CodlMsHui +£'’*(D)) = Co)! + £'’*(D)) by Lemma FTlT bL one 
obtains 

-^$(«.) > m^^fsWh - c\\y^fs\\lt^ - cwv^fsh^ > lps\\y^M\h - C', 


(7.12) 












14 


DANIEL LOIBL, DANIEL MATTHES, AND JONATHAN ZINSL 


which is bounded from below, using Young’s inequality in the last step. We are now concerned with the 
passage to the limit inferior as s 0. In view of the flow interchange lemma fTheorem l2.2l) . we can assume 
that at least for small s > 0 the quantity --^^{us) is bounded from above. By (I7.12|) . we infer that 
I|V^/s||l 2 is bounded in s. We conclude that (on a suitable subsequence) f{us) /(u?) in as well as 

f{us) —>■ /(u”) in by Rellich’s theorem, as s \ 0. By continuity of /, one also has Ug —>■ u" pointwise 

almost everywhere on fl. By a straightforward application of Vitali’s theorem and Lemma |7.1I cl. one gets 
that f'{us)us —>■ /'(u”)m” in Hence, weak-strong convergence leads to 

lim f [Af{us)Vfs-^(f + Afsf{ug)usA(p] dx = [ [A/(-«”)V/«) • Vv? + dx 

= [ A/>«,(p)dx. 

Jn 

We thus arrive at 

2 )'^($«))> [ Mf{K,^)dx. 

Jn 

The flow interchange lemma fTheorem l2.2p then yields for ip and —ip, respectively: 




- -w; 


O [ K ^log< ^dx + P f <logu"da; 

^ Jn Jn 

n—1 


< / u'!^ipdx— I u" ^ipdx + T f Aff , ip) dx 

Iq J q J q 


< Sw2«,<-1)+/3 [ <-ilog<-Mx-/3 [ <log<da:. 

^ Jn Jn 

Proceeding similarly as for Lemma ISTTl we introduce a nonnegative test function rj G C))°((0, oo)) fl C(R+), 
multiply the chain of inequalites above with 77 (nr) and sum over n G N to eventually get 

f r]r{t) - T]r{t + t) 


< 


- lhllcoKT$(Mo) +/3 

Jo Jn 

Vrit) -r]r{t + t) 


-Ur log(uT-) da; dt 


Jo Jn ^ 

< \\r]\\coKT^{uQ) - P 


Uripdxdt-\- / / rjrJJ'f {Ur, ip) dx dt 

Jo Jn 

' f Vr{t) - r]r{t -G t) 


(7.13) 


0 Jn 


-Ur log(MT) da: dt. 


To arrive at (TTtI) for arbitrary test functions 77 G C))°((0, 00 )) fl C(]R+), we decompose r] = r]+ — rj- into its 
positive and negative parts and use (17.131) for rj^ and rj-, respectively (recall that jTyj = 77 + -f 77 _). □ 


As in the previous section, the proof of Theorem II.21 is completed by passing to the continuous-time limit 
r \ 0. 

Proposition 7.6 (Continuous-time limit). Let T > 0, a vanishing sequence of step sizes rj, 0 (fc —> 00) 
and the associated family of discrete solutions {uri^)keN given. Then, the following statements hold with 
p > 1 from Proposition \ 7. 5\ 

There exists a (non-relabelled) subsequence and a map u G C'^/^([0,r]; (.i^(O), W 2 )) D L°°{[0,T]; Lp{Q)) 
such that for fc —>■ 00; 

(a) Ur,, —>• u uniformly w.r.t. t G [0,r] in {3^{Tt),'W 2 ) ■ 

(b) Ur,. u in LP(\d,TV LP(LI)) and almost everywhere on [0,Tl x fl. 

(c) f{ur,)^f{u) tnL\[0,T];H^{n)). 

(d) f{ur,)^fiu) ^nL^i[0,T]■,H^ifl)). 

The limit u is a weak solution to dni) in the sense stated in Theorem \1.4\ and for almost every t G [0,T], 
one has 


lim <t>{ur,,{t,-)) = <P{u{t,-)). 

k—yoo 

Proof. Clearly, /(ur*,) ^ u in L^([0, T]; i7^(H)) for some limit v by Alaoglu’s theorem and estimate (17.81) . We 
show that f{ur) —>■ f{u) in L^([0, T]; L^(H)) for u = f~^{v). Then, by a standard interpolation inequality. 
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f{ur) f{u) in L^([0, T]; iJ^(f2)) follows. We seek to apply Theorem 12.31 and let X := LP{Q) (with p > 1 
from Proposition [7^, 


A[p) := 


WfipWm. 

+ 00 , 


if/(p) 6^1(0), 
otherwise, 


and g as in the proof of ProDOsition l6.2l Obviously, •)) dt is A:-uniformly bounded thanks to (17.81) . 

For the application of Theorem 12.31 it remains to verify that A has relatively compact sublevels in X. Let a 
sequence [pk)km in ^ with A{pk) < C for all fc G N and some C > 0 be given. Clearly, \\f{pk)\\H^ < Vc 
for all k; hence, by the Rellich-Kondrachov compactness theorem, f{pk) f{p) in U{Vl) for sufficiently 
large r G (l) ns well as pointwise a.e. on fl, for some limit f{p) G L''{n), extracting a suitable 

subsequence. By continuity, pk ^ p pointwise a.e. on 17. Without loss of generality, we may also assume 
that pk —>■ p weakly in X, by Alaoglu’s theorem and dLH) — which also yields L^’-uniform integrability of 
(pfe)fcGN- Convergence of pk to p in X then follows by Vitali’s convergence theorem. 

The application of Theorem 12.31 yields the existence of a limit map u such that (on a subsequence) 
—>■ in X, in measure w.r.t. t G [0,T]. As in the proof of Proposition 16.21 we deduce that 

Mrs, —>■ M in L^’([0, T]; LP(17)) and pointwise almost everywhere on [0,T] x 17, possibly extracting further 
subsequences. Hence, by continuity, f(ur^) —>■ f(u) a.e. on [0, T] x 17. Since is fc-uniformly bounded in 

L°“([0,T];(17)) thanks to Sobolev’s inequality and (17.8p . the family (/(uT-fc))fcGN is uniformly integrable 
in T^([0, T]; T^(17)). Vitali’s theorem yields the claimed convergence in L^([0, T]; L^(17)). 

The limit u is a weak solution to (EB because of the following: First, since — ^ < zlogz < z^ for 
all z > 0 and thanks to the uniform estimate in L°“([0,T];LP(17)), Ur^ log(MTfc) is fc-uniformly bounded in 
i^([0, T]; T^(17)). Hence, letting k ^ oo and /3 \ 0 in (17.7|) yields 


lim 

k—¥oo 



■nrkit + rk) 
Tk 


Ut^ Lp dx dt 



VTk-^fiurk,p) dx dt 


= 0 . 


We obtain the time-continuous weak formulation since 

poo p poo p 

lim / / 77T-fcA//(urfc, p) dxdt = / / pNf{u,ip) dx dt, 

dn do Jo. 

thanks to the weak convergence Af^Ur^) Af{u) and the strong convergence yf^Ur^) —?• V/(u), both 
in T^([0, T]; L^(17)). Especially, Ur,,f'{ur,,) uf'{u) in T^([0, T]; L^(17)) due to Lemma ITT] and Vitali’s 
theorem: 

f [ {uTj'{UrJ)'^dxdt<Ci [ [ {f{UrJ + l)'^dxdt<2Ci [ [ {f{Urj'^ + 1) dx dt, 

Jo Jo Jo Jo Jo Jo 

and (/(uTfc))feGN is uniformly integrable in L^([0, T]; L^(17)). □ 


Appendix A. 

Lemma A.l (Binomial formula for traceless symmetric matrices). Let A G be symmetric and traceless. 
Then, for all v G M'’*, the following estimate holds: 

llAf+ 2u->0. (A.l) 

d 

Proof. Without loss of generality, |u| = 1. Using the orthogonal decomposition A = SDS"^ for orthogonal S 
and diagonal D, (lA.ll) is — by the transformation w := S'^v (note that |w| = 1) — equivalent to 


XI + 2 ^ A, 


d- 1 


w. 


2=1 


2 = 1 


d 


> 0 , 
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where Ai,..., are the (real) eigenvalues of A. Without restriction, assume that Ai is the smallest eigenvalue 

d d 

of A. Since A is traceless, we have X! Since Aiwf > Ai, we obtain 

2=1 2=1 


d d j ^ d / d \ d 

>^^Wi-\ -+ 1 -251^*-I- T~ 


2 = 1 


2 = 1 


2=2 


^2 = 2 


2=2 


= E A.- 


2=2 


E A.- 


> 0 . 


□ 
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